Shape transformations of a compartmentalized fluid surface 
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A surface model on compartmentalized spheres is studied by using the Monte Carlo simulation 
technique with dynamical triangulations. We found that the model exhibits a variety of phases: 
the spherical phase, the tubular phase, the planar phase, the wormlike planar phase, the wormlike 
long phase, the wormlike short phase, and the collapsed phase. It is also shown that almost all 
phases are separated from their neighboring phases by first-order transitions. Mechanical strength 
of the surface is given only by elastic skeletons, which are the compartment boundaries, and vertices 
diffuse freely inside the compartments. We confirm that the cytoskeletal structure and the lateral 
diffusion of vertices are an origin of such a variety of phases. 

PACS numbers: 64.60.-i, 68.60.-p, 87.16.Dg 



I. INTRODUCTION 

Biological membranes and synthetic polymer- 
membranes show a variety of shapes such as spherical, 
tubular, discoid, cylindrical, and many others including 
starfish The shape of membranes is partly under- 
stood numerically with a surface model called a minimal 
model i^] and also with the area difference bilayer model 
[3, External forces such as gravity and flow fields 
make the surface shape change US 0| ■ 

The surface shape can also be influenced by thermal 
fluctuations [1, Therefore, the membrane shape 

should be understood as an equilibrium statistical me- 
chanical phenomenon, although the shape of membranes 
seems to have a non-equilibrium nature. We should re- 
mind ourselves that the shape transformation and the 
surface fluctuation are two different phenomena, where 
the surface fluctuation phenomena have been extensively 
studied statistical mechanically [ll|, El [H, 13, [H, . 

The surface fluctuation transition is accompanied by 
the collapsing transition in artificial membranes and 
in the surface model [H, ^ of Helfrich Polyakov 
[2l| , and Kleinert [l^l ■ This also indicates that the shape 
of membranes should be understood within the context 
of the theory of phase transitions. In fact, shape trans- 
formations, such as the prolate-oblate transition driven 
by the thermal fluctuation, were experimentally observed 
[31 • Current understanding of the effects of thermal 
fluctuations on the conformation and the elastic prop- 
erties of membranes are reviewed in 10]. In [23|, it was 
shown that the shape of a compartmentalized fluid sur- 
face model changes due to thermal fluctuations. The 
results suggest that possible origins for the variety of 
membrane shapes are the cytoskeletal structure and the 
fluidity of lipids in membranes. Moreover, it was also 
suggested in [2^ that the large variety of shapes can be 
understood in the framework of a surface model which 
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has a cytoskeleton. The cytoskeletal structure has been 
considered as a key notion for understanding physics of 
membranes [M Hi [H ■ 

In this article, in order to make this observation more 
convincing we study another compartmentalized fluid 
surface model, which is almost identical to the model 
in It is remarkable that a small change in the 

model makes a large difference in the multitude of surface 
shapes. The only difference between the model in this ar- 
ticle and that of Ref . [23j is in the junction elasticity; rigid 
plates are assumed as the junctions in [11] while neither 
two-dimensional elasticity nor rigid plate is assumed in 
the model of this article. Both of the compartmental- 
ized models are inhomogeneous because of the cytoskele- 
tal structures; the surface strength on the compartment 
boundary is different from that inside the compartments, 
and moreover the diffusion of vertices is confined only in 
the compartments. 



II. MODEL 

The compartmentalized structure is a sublattice on a 
triangulated surface, which is constructed from the icosa- 
hedron. By dividing the edges of the icosahedron into l 
pieces, we have a triangulated lattice of size N = 10£^ -1-2, 
which is the total number of vertices. Then, we have a 
sublattice of size Ns = 30m£ in the N = 10^^-1-2 lattice 
if m divides t. The vertices in the sublattice include the 
junctions of the compartments on the N — 10^^-1-2 lat- 
tice, and the total number of junctions Nj is given by 
Nj = 10m^-|-2. The total number of links between the 
junctions is 3A^j— 6, and each link contains £/m vertices. 
Thus we have Ns = 30m£. The compartment size can 
be characterized by ti = X^iffi"'*"^ which is the total 
number of vertices in a compartment. 

Figure[TJa) shows the starting configuration for Monte 
Carlo simulations. The size of surface is characterized 
by two integers (£, m) = (16, 2), and the size is given by 
(A^,7Vs,7Vj) = (2562,960,42), and n = 21. 

The model is defined by the Gaussian bond potential 
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FIG. 1: (Color online) (a) Starting configuration of surfaces 
of size {N, Ns,Nj) = (2562, 960, 42); thick lines denote a sub- 
lattice composed of the linear chains and the junctions, (b) 
tangent vectors ti and t2 at a vertex with coordination num- 
ber q = 6 that give rise to a contribution 1 — ti • t2 to the 
bending energy with the weight of 1, and (c) tangent vectors 
ti, t2 and ts at a vertex with coordination number q = 5 that 
contribute 1 — [ti ■ (t2-|-t3)]/2 to the bending energy with the 
weight of 1/2. 



Si and the one-dimensional bending energy S2 , which are 
respectively defined by 
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where Xi is the three-dimensional position of the vertex 
i and is a unit tangent vector of the bond i. 

is the sum over all bonds {ij) on the lattice, and J2{ij) 
in S2 is the sum over all nearest neighbor bonds (ij) on 
the sublattice. 

Tangent vectors at the junctions of coordination num- 
bers q — 6 and q — 5 are shown in FigslH^b) andfljc). The 
tangent vectors ti and t2 at a junction of coordination 
number q — 6 give rise to a contribution 1 — ti • t2 to ^2 
with the weight of 1. The remaining two inner-products 
of tangent vectors are defined just like 1— ti • t2 at the 
q = 6 vertices. On the contrary, the tangent vectors ti, 
t2 and ta at a junction with coordination number q = 5 
contribute 1— [ti • (t2-f ts)] /2 to the bending energy with 
the weight of 1/2. The remaining four inner-products of 
tangent vectors are defined just like 1 — [ti • (t2-|-t3)] /2 
at the q = 5 vertices; this is the reason for the weight 1/2 
of the bending energy 5*2 at the q = 5 junctions. Conie- 
quently, the definition of the bending energy at the q — 6 
junctions is almost identical to that at the q — 5 junctions, 
whose total number is only 12. 

The partition function Z of the model is given by 



// « 
l[dX,exp[-S{X,T)], 

1 i=l 

S{X,T) = Si + bS2, 
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where S{X, T) is the Hamiltonian, and b[kT] is the 
bending rigidity, which is a microscopic quantity and 
therefore, it is not always identical to the macroscopic 
bending rigidity. denotes all possible triangulations 
which keep the compartment boundary (= the sublat- 
tice bonds) unchanged. /' Oi^i d^Xi denotes the multiple 
three-dimensional integrations under the condition that 
the center of mass of the surface is fixed. 



III. MONTE CARLO TECHNIQUE 

The integrations of the dynamical variables X and 
T are performed by the Monte Carlo simulation tech- 
nique [13, m, [s^, [3l|. The three dimensional ran- 
dom shift SX of X generates a new position X' = X+X , 
which is accepted with the probability Min[l, exp(— A5)], 
where AS = S{new) — S{old). The vertices can be clas- 
sified into three groups; the vertices inside the compart- 
ment, the vertices on the compartment boundary, and 
the vertices at the junctions; the final two groups of ver- 
tices are those constructing the sublattice. The point 
dX is randomly chosen in a sphere and, the radius of 
the sphere is fixed at the beginning of the simulations 
so that the acceptance rate is equal to about 50% in 
each group of vertices. The radius assumed for one 
group of vertices is not always identical to those for the 
other groups of vertices. The summation over T is per- 
formed by the standard bond flip technique [1^, [sO, |3l| 
and, therefore the acceptance rate for the flip is not 
fixed a priori but found to vary approximate bracket 
70% ~ 75%, which is slightly dependent on b. We as- 
sume the surfaces of size {N,Ns,Nj) = (5762,2160,92) 
and {N,Ns,Nj) = (10242,3840,162), which correspond 
to integers {£, to) = (24, 3) and {£, m) = (32, 4). The total 
number of MCS (Monte Carlo sweep) after the thermal- 
izaion is about 1x10* '-^ 1.5x10* for the A^ = 5762 surface 
and 1.3x10* - 2x10* for the 10242 surface. The ther- 
malization process comprises about 1 x 10* MCS, which 
is sufficiently large, in almost all cases. 



IV. RESULTS OF SIMULATION 
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FIG. 2: The mean square size versus b at (a) small b region 
b<20 and at (b) whole region 6<240. The compartment size 
is given by n = 21, which is the total number of vertices in a 
compartment. 

The shape of surfaces can be reflected in the mean 
square size X^, which is defined by 
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where X is the center of the surface. Figures [HJa) and 
[DJb) show X'^ versus b obtained at relatively small b 
region and at whole b region, respectively. The verti- 
cal dashed lines denote phase boundaries where X^ dis- 
continuously changes. We have at least seven phases 
in the region < 6 < 240 on the surface of size 
{N,Ns,Nj) = (10242,3840,162). We cah the phases 
as the collapsed, the wormlike short, the wormlike, the 
wormlike planar, the planar, the tubular, and the spher- 
ical. Almost all of two neighboring phases, except the 
wormlikc planar phase and the planar phase, seem to be 
connected by a first-order transition, because X^ discon- 
tinuously changes at the phase boundaries. 

X^ in the wormlike planar phase of the 10242 sur- 
face is wildly fluctuating on the b axis; X^ at b=15 and 
those at the phase boundary close to the planar phase are 
different from the remaining X^ in the wormlike planar 
phase. X^ at b—13 of the N = 5762 surface also seems 
to be an anomalous value. The surface shape at 6=15 of 
the iV = 10242 surface and that at &=13 of the iV = 5762 
surface are wormlike, and then we understand that the 
configuration was trapped in the potential minimum cor- 
responding to the wormlike long phase in the simulations. 
The potential barriers separating the phases seem to be 
low because the surface size is not sufficiently large, and 
for this reason such anomalous behavior of X'^ can be 
seen in the wormlike planar phase. We must emphasize 
that the anomalous behavior of X"^ does not imply that 
the model is ill-defined. In fact, the Hamiltonian such as 
the bending energy ^2 is not unstable and has the unique 
value corresponding to the given value of b even in the 
wormlike planar phase as we will see later in this paper. 

Snapshots of surfaces and their sections are shown in 
FigsI2Ja)-[3fe), which were respectively obtained in the 
collapsed phase, the wormlike long phase, the planar 
phase, the tubular phase, and the spherical phase. The 
surfaces and the surface sections were shown in the same 
scale. The self-avoiding property [H, [sj is not as- 
sumed in our model and, therefore the phase structure in 
the small b region seems phantom. However, as we see in 
the snapshots, the phase structure seems realistic in the 
large b region. 

In order see the difference between the wormlike short 
phase, the wormlike long phase, and the wormlike pla- 
nar phase, we show snapshots of skeletons in FigsHJa)- 
SJd), which were obtained in the wormlike short phase 
(FiglHa)), the wormhke long phase (FigH^b)), and the 
wormlike planar phase (FigsUJc) andUI^d)). All figures 
were drawn in the same scale. We understand from the 
snapshots in FigsHJc) andjU^d) that one part of the sur- 
face is wormlike and the remaining part is planar in the 
wormlike planar surfaces. We see from FigsHI^c) and|D^d) 
that the size of planar part varies depending on b in the 
wormlike planar phase; we see the inflated parts are pla- 
nar from their surface sections. It is easy to understand 
that X"^ is strongly dependent on the size of the planar 
part. X'^ also depends on the number of planar parts; 
we see two planar parts at the two ends of the surface 



in Figllljd). For this reason, X"^ wildly fluctuates at the 
phase boundary (6 ~ 35) between the wormlike planar 
phase and the planar phase as mentioned above. 

The surface shape in the wormlike short phase is also 
wormlike as we see in Fig|4l[a), however, the thickness or 
equivalently the longitudinal length of surface is slightly 
different from those of the surfaces in the wormlike long 
phase. This difference is reflected in X"^, and conse- 
quently the wormlike short phase is separated from the 
wormlike long phase by the flrst-order transition. 

The one-dimensional bending energy S2/N'g versus b 
is shown in FigsOJa) and [5{b), where N'g is the total 
number of vertices where ^2 is defined. The junctions 
of coordination number q = 6 (9 = 5) are counted 3 (2.5) 
times in N'g because of the definition of S'2 and, therefore 
N'g is given by Ng = Ns + 2Nj — 6, which is also written 
by N'g = 10£^-60m^ + 2. The vertical dashed lines in 
the figures denote the phase boundaries. We find a dis- 
continuous change in S2/N'g at the boundaries between 
the collapsed phase and the wormlike short phase and at 
the boundary between the wormlike short phase and the 
wormlike long phase. No discontinuous change can be 
seen in S2/Ng at any other phase boundaries. 

We see no wild fluctuation of S2/N'g in the wormlike 
planar phase in FigslSfa) and[5{b). S2/N'g smoothly vary 
even at 6=15 and at 6~35, where X'^ wildly fluctuates. 

The two-dimensional bending energy is defined by 

^3 = E (1 - • ' (4) 

(y) 

where n.^ is a unit normal vector of the triangle i. The 
surface fluctuation can be reflected in S^, which is not 
included in the Hamiltonian. FiguresjUJa) and[SKb) show 
Ss/Nb versus 6, where Nb = 5N — 6 is the total number 
of bonds. Discontinuous changes can be seen in S^/Nb 
at the boundary between the wormlike short phase and 
the wormlike long phase, at the boundary between the 
wormlike long phase and the wormlike planar phase, at 
the boundary between the planar phase and the tubular 
phase, and at the boundary between the tubular phase 
and the spherical phase. The discontinuous changes in 
Ss/Nb are consistent to those in X^ in Figsl^a) and 
[2ljb). Note also that anomalous spikes of Ss/Nb at 6= 15 
and at 6 = 35 correspond to the anomalous value or the 
wild fluctuations of X^ mentioned above. 

We have seen that at least one physical quantity dis- 
continuously changes at the phase boundaries except the 
boundary between the planar phase and the wormlike 
planar phase. At this boundary we see that no physical 
quantity discontinuously changes although AT^ anoma- 
lously fluctuated, which was seen in Fig. [Ijb). Then, the 
discontinuous nature of the transition at this boundary 
is not confirmed from the numerical data in this paper. 
Therefore, we consider that almost all phases, except the 
planar and the wormlike planar phases, are separated 
from their neighboring phases by first-order transitions. 

The Gaussian bond potential Si/N is expected to be 
S'l/A^ ~ 3/2 because of the scale invariance of the par- 





(a) 6 = 6 (b)&=10 (c)6=14 (d) 6 = 35 

FIG. 4: (Color online) The snapshots of skeletons obtained at (a) 6 = 6 (wormlike short phase), (b) 6 = 10 (wormlike long 
phase), (c) 6=14 (wormlike planar phase), and (d) 6 = 35 (wormlike planar phase). The surface size is given by {N, Ns , Nj) — 
(10242,3840,162). 
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FIG. 5: The bending energy S2/N'g versus 6 at (a) small 6 
region 6<20 and at (b) whole region 6<240. 
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FIG. 6: The two-dimensional bending energy Ss/Nb versus 6 
at (a) small 6 region 6 < 20 and at (b) whole region 6 < 240. 



tition function. Figures [7fa) and [7{b) show that the 
expected relation is almost satisfied. We find that the 
relation is satisfied in the region of low bending rigidity, 
where the surfaces are almost collapsing. On the con- 
trary, we find that the relation is not exactly satisfied in 



the region of high bending rigidity, where the surfaces are 
inflated, although the deviation is very small compared 
to the value itself. 
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FIG. 7: The Gaussian bond potential Si/N versus b at (a) 
small b region 6 < 20 and at (b) whole region b < 240. 



SUMMARY AND CONCLUSION 



the spherical phase, the tubular phase, the planar phase, 
the wormlike planar phase, the wormlike long phase, the 
wormlike short phase, and the collapsed phase. Almost 
all two neighboring phases are connected by first-order 
transitions. The spherical phase and the tubular phase 
are connected by a first-order transition, which is quite 
similar to the prolate-oblate transition. Our results in- 
dicate that the variety of membrane shapes and their 
transformations can be understood in the inhomogencous 
model, which is characterized by compartmentalization 
of fluidity of vertices and the cytoskeletal structure con- 
structed on the conventional homogeneous surface model. 

It is interesting to study the model by including the 
two-dimensional bending energy in the Hamiltonian. The 
phase structure in the thermodynamic limit and the de- 
pendence of the phase structure on the compartment size 
still remains to be clarified. 



To summarize the results, we have investigated a com- 
partmentalized fluid surface model by using the canonical 
MC simulation technique and found a variety of phases; 
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